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A STUDY OF INVISCID FLOW ABOUT AIRFOILS AT HIGH SUPERSONIC SPEEDS ! 


By A. J. Есаввв, Jr., CLARENCE A. SvvERTSON, and Samus. Kraus 


SUMMARY 


Steady flow about curved airfoils at high supersonic speeds is 
investigated analytically. With the assumption that air behaves 
as an ideal diatomic gas, it is found that the shock-expansion 
method may be used to predict the flow about curved airfoils up to 
arbitrarily high Mach numbers, provided the flow deflection 
angles are not too close to those corresponding to shock detach- 
ment. This result applies not only to the determination of the 
surface pressure distribution, but also to the determination of the 
whole flow field about an airfoil. Verification of this observa- 
tion is obtained with the aid of the method of characteristics by 
extensive calculations of the pressure gradient and shock-wave 
curvature at the leading edge, and by calculations of the pressure 
distribution on a 10-percent-thick biconvex airfoil at 0° angle of 
attack. 

An approximation to the shock-expansion method for thin 
airfoils at high Mach numbers is also investigated and is found 
to yield pressures in error by less than 10 percent at Mach 
numbers above 3 and flow deflection angles up to 259. This 
slender-airfoil method is relatively simple in form and thus voy 
prove useful for some engineering purposes. 

Effects of caloric imperfections of air manifest n. ied 
flow fields at hagh Mach numbers are investigated, particular 
attention being given to the reduction of the ratio of specific 
heats. So long as this ratio does not decrease appreciably below 
1.8, 18 indicated that the shock-expansion method, generalized 
to include the effects of these imperfections, should be substan- 
tially as accurate as for ideal-gas flows. This observation is veri- 
fied with the aid of a generalized shock-expansion method and a 
generalized method of characteristics employed in forms appli- 
cable for local air temperatures up to about 5,000° Rankine. 

The slender-airfoil method is modified to employ an average 
value of the ratio of specific heats for a particular flow field. 
This simplified. method has essentially the same accuracy for 
imperfect-gas flows as its counterpart has for ideal-gas flows. 

An approximate flow analysis is made at extremely high Mach 
numbers where it is indicated that the ratio of specific heats may 
approach close to 1. In this case, it is found that the shock- 
expansion method may be in considerable error; however. the 
Busemann method for the limit of infinite free-stream Mach 
number and specific-heat ratio of 1 appears to apply with 
reasonable accuracy. ‘ 


INTRODUCTION 


Small-disturbance potential-flow theories have been em- 
ployed widely, and for the most part successfully, for pre- 


dicting the pressures (and velocities) at the surface ofan 
airfoil in steady motion at low supersonic speeds. Thus the 
linear theory of Ackeret (ref. 1) has proven particularly use- 
ful in studying the flow about relatively thin, sharp-nosed 
airfoils at small angles of attack, while the second-order 
theory of Busemann (ref. 2) has found application when 
thicker airfoils at larger angles of attack were under consid- 
eration. At highfree-stream Mach numbers the range of appli- 
cability of any potential theory is seriously limited, however, 
due to the production of strong shocks by even the relatively 
small flow deflections caused by thin airfoils. The assump- 
tion of potential flow is invalidated, of course, by the pro- 
nounced entropy rises occurring through these shocks. 

This limitation on potential theories was early recognized 
and led to the adoption (see ref. 3) of what is now commonly 
called the shock-expansion method. The latter method de- 
rives its advantage over potential theories, principally, by 
accounting for the entropy rise through the oblique shock 
emanating from the leading edge of a sharp-nosed airfoil. 
Consequently, so long as the disturbed air behaves essen- 
tially like an ideal gas, and so long as entropy gradients 
normal to the streamlines (due to curvature of the surface) 
do not significantly influence flow at the surface, the shock- 
expansion theory should predict the pressures at the surface 
of an airfoil with good accuracy—it is tacitly assumed, of 
course, that the flow velocity is everywhere supersonic, and 
that the Reynolds number of the flow is sufficiently large to 
minimize viscous effects on surface pressures. 

The departure of the behavior of air from that of an ideal 
gas at the temperatures encountered in flight at high super- 
sonic speeds has been the subject of some investigation in 
the case of flows through oblique shock waves. In reference 
4, the effects of thermal and caloric imperfections on the 
pressure rise across an oblique shock wave were investigated 
at sea-level Mach numbers of 10 and 20, and it was found 
that these effects decreased the rise by less than 5 percent 
for maximum temperatures up to 3,000° R. (corresponding 
to flow deflection angles up to 24°). This decrease was found 
to be due almost entirely to caloric imperfections, or changes 
in vibrational heat capacities of the air passing through the 
shock wave. The changes in temperature and density of the 
air passing through the wave were affected to a considerably 
greater extent. Subsequently, an investigation was carried 
out by Ivey and Cline up to Mach numbers as high as 100 
(ref. 5), using the results for normal shock waves obtained 
by Bethe and Teller considering effects of dissociation (ref. 


1 Supersedes МАСА TN 2646, “Inviscid Flow About А1015 at High Supersonic Speeds” by A.J. Eggers, Jr., and Olarenco A. Syvortson, 1952, and NAOA TN 2729 “Ап Analysis of 


Supersonic Flow in the Region of the Leading Edge of Curved Atrfoils, Including Charts for Determining 


Surtace-Prossure Gradient and Shock-Wave Curvature” by Samuel Kraus, 1952, 
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6). As would be expected, the pressures were found to be 
affected to a somewhat greater extent at the higher 1 Mach 
numbers. 

The extent to which flow in the region of the leading edge 
of an airfoil departs from the simple Prandtl-Meyer type 


has also been investigated at high supersonic airspeeds. If 


the surface is curved, for example, to give an expanding 
flow downstream of the leading edge, expansion waves from 
the surface will interact with the nose shock wave, thereby 
curving it and yielding a nonisentropic flow field. This 
flow field may be characterized not only by disturbances 
emanating from the surface but also by disturbances reflect- 
ing to some extent from the shock wave back toward the 
surface. The manner in which these phenomena dictate 
shock-wave curvature and surface pressure gradient in 
ideal-gas flows at the leading edge has been treated by 
Crocco (ref. 7) and more recently by Schaefer (ref. 8), Munk 
and Prim (ref. 9), Thomas (ref. 10), and others. Іп the 
cases considered by Munk and Prim, it was found that sur- 
face pressure gradients were less (in absolute value) than 
those obtained assuming Prandtl-Meyer flow at the higher 
Mach numbers (i. e., Mach numbers greater than about 3) 
although, generally, by no more than about 10 percent. 
Since curved airfoils are likely to be of fundamental interest 
at high flight speeds (see, e. g., ref. 11), the effects of reflected 
disturbances would appear to merit further investigation, 
both at the leading edge and as regards their influence on 
the whole flow field. In addition, it would appear desirable 
to consider effects of gaseous imperfections through the field. 

Such an investigation has therefore been undertaken in 
the present report, using the method of characteristics to 
obtain accurately flow fields and as a basis for obtaining the 
more approximate methods of analysis. The method is 
employed in a generalized form which allows caloric imper- 
fections, as well as entropy gradients, in the flow to be con- 
sidered at temperatures up to the order of 5,000° R.—thermal 
imperfections are neglected (see ref. 4). A 10-percent-thick 
biconvex airfoil is treated at Mach numbers from 3.5 to 
infinity, and the results are compared with the predictions 
of the shock-expansion method, including & simplified form 
of the method applicable to slender airfoils at high Mach 
numbers and a generalized form of the method including 
effects of caloric imperfections. In addition, flow’ in the 
region of the leading edge of curved airfoils is considered in 
some detail. Values of the surface pressure gradient and 
shock-wave curvature are presented for a wide range of 
Mach numbers and flow deflection angles. | 


SYMBOLS 

a local speed of sound 

с chord 

C, С; characteristic coordinates (Сі positively inclined 
and О, negatively inclined with respect to the 
local velocity vector) ` 

Ca section drag coefficient 

Ci section lift coefficient 

Cm section moment coefficient (moment taken about 


leading edge) 
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С, pressure coefficient, ага 
ба specific heat at constant pressure 
б; specific heat at constant volume 
К curvature 
М Mach number (ratio of local velocity to local speed 
of sound) 
Р pressure ratio, Р. 
Po 
p -  gtatic pressure 
q dynamic pressure 
R gas constant 


8,7 rectangular coordinates (in streamline direction 
and normal to streamline direction, respec- 


tively) 

Г temperature, °R 

і time 

V resultant velocity 

W distance measured from leading edge along airfoil 
surface 

x, Y rectangular coordinates 

а ` angle of attack, radians unless otherwise specified 

Mach angle, arc sin (sr , radians 

Y ratio of specific heats, z 
(Average value of y is Ya.) 

8 flow deflection angle, radians unless otherwise 
specified | 

t angle between shock wave and flow direction just 
downstream of shock wave, radians 

0 molecular vibrational energy constant, °R. (5,5009 
R. for air) 

K ratio of shock-wave curvature to that given by the 


shock-expansion method 


p mass density 
т ` ghock-wave angle, radians 
y ratio of surface pressure gradient to that given by 
the shock-expansion method 
د‎ ray angle for Prandtl-Meyer flow, radians 
SUBSCRIPTS 
О `-  free-stream conditions 
A,B,C,D conditions at different points in flow field 
Ü ideal-gas quantities 
N ‘conditions at the leading edge immediately down- 
` stream of the shock wave 
S conditions on streamline 
10 conditions along airfoil surface 
g^ conditions along shock wave 
SUPERSCRIPT 


- - vector quantities 
DEVELOPMENT OF METHODS OF ANALYSIS 
GENERAL METHODS 


Method of characteristics. —Two-dimensional rotational 
sup ersonic flows have been treated by numerous authors with 
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the aid of the method of characteristics, and various adapta- 
tions of the method have been found which are especially 
suited for studying particular types of such flows. In the 
case of steady flows in which atmospheric air does not behave 
as an ideal diatomic gas, a very familiar and simple form of 
the compatibility equations may be employed. То illustrate, 

consider the Euler equation 


p dV... grad 2 (1) 


the continuity equation 


div (pV)=0 (2) 


and the equation for the speed of sound (evaluated at con- 
stant entropy) 


d 


| )3( لک 


Rewriting equations (1) and (2) in the form of partial dif- 
ferential equations and transforming the resulting expressions 
to the characteristic, or Cj, С, coordinate system, there is 
obtained, upon combination with equation (3), the following 
relations for steady flow: 


cot vi Op m 
ТӨЛІ кіні (4) 
and | 
cot В Ор 
pV? sed 35.) (56.96. (5). 


A simple addition or subtraction of equations (4) and (5) 
then yields the compatibility equations (see, e. g., ref. 12) 


Op | ò — 
aa? tan B so (6) 
and ` 5 5 
D E 
5б ^V tan B 20; (7) 


Now, іп reference 4 both caloric and thermal imperfections 
of air were considered, and it was found that the latter 
imperfections 2 have a negligible effect on shock processes in 
atmospheric air. It may easily be shown that this conclusion 
also applies to expansion processes and, for this reason, 
caloric imperfections only are considered in detail in the 
present paper. These imperfections become significant in 
air at temperatures greater than about 800° R. and first 
manifest themselves аз changes in the vibrational heat 
capacities with temperature. Thus, the specific heats c, and 
с, and their ratio y for the gas also change., The equation of 
state remains, however, 


p=phT (8) 


and the specific heats are still related to the gas constant by - 


the expression 
€, —C,— № (9) 


2 Thermal] imperfections usually appear In thc form of intermolecular forces and molecular- 
size effects and may be accounted for with additional terms in the equation of state. 
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Furthermore, it readily follows from the differential energy 
equation and these expressions that the speed of sound is 
given by the simple relation 


G—yRT | (10) 
Combining equations (8) and (10) and noting that 
sin 8—2a/V there is then obtained 
_ Үр? 


Hence, on combining this equation with equations (6) and 
(7), it is apparent that the familiar compatibility equations 


d0, sin 28 20, (12) 
and 
a Op. 2үр 29 


also hold for the more general type of flow under considera- 
tion. ‘These equations are basic, of course, to two-dimen- 
sional characteristics theory, and, as will be shown later, 
form а convenient starting point for developing simpler 
theories of two-dimensional supersonic flow. 

In order to apply equations (12) and (13), it is evident 
that the manner in which y and В.ог M are connected to 2 
or 6 must be known, Relations implicitly connecting these 
variables at temperatures up to the order of 5000° R may be 
readily obtained from the results of reference 4 by simply 
eliminating the terms therein accounting for thermal imper- 
fections. Thus we have as a function of the local static 
temperature and free-stream conditions 


D 1) (2) әгі etm . (14) 


1+(y:—1) Ө) (eT —1)? 1)? 


GP es) amaze) | 
(15) 


For isentropic flow along a streamline, the pressure is related 
to the temperature by the expression 


Ү— 4 


апа 


p, AQ) 
| p», A(T) ы 
where 
EN dies e?/T — 1 E 
ا‎ 0 


If there is a shock wave in the flow,? in particular, a nose or 
leading-edge shock, then the follówing additional relations 
obtained with equations (8), (10), and (15) and the conditions 


3 If there are no shock waves, then the subscript о in-equation (10) can, of course, be replaced 
with the subscript 0. 
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for continuity of flow and conservation of momentum along 
а streamline through the shock are also required: 


(MA ам 


y | а+=мэ—р@ Md | +4 T азу 


ы t Ie TM os MA1 
sin! o= (ҰЯ Гасан (19) 
апа 
tan a С E (2 0) 
(ps[po) — 1 


By use of the local static temperature as & parameter, the 
term 2yp/sin 28 in equations (12) and (13) may now be 
evaluated with equations (14) through (17). Equations (18) 
through (20) define the initial conditions downstream of a 
leading edge or other shock wave in the flow field. Thus, 
equations (12) through (20) provide all the information 
necessary to calculate the flow about an airfoil by, means of 
the method of characteristics. As described in detail in 
Appendix A, the calculations are of three general types: 
(1) determination of conditions at a point in the flow field 
between the shock and the surface; (2) determination of 
conditions at a point on the surface; and (3) determination 
of conditions at a point just downstream of the shock. Case 
(1) entails the use of both compatibility equations, while 
case (2) entails the use of the compatibility equation for a 
second-family characteristic line in combination with the 
equation of the airfoil surface, and case (3) involves the 
compatibility equation for a first-family line їп combination 
with the oblique-shock equations.. With the aid of the three 
general types of calculations, the entire flow field about an 
airfoil сап be built up numerically using a computing pro- 
cedure working from the leading edge downstream. Іп cases 
where changes in the vibrational heat capacities with tem- 
perature are neglected, the calculations are, of course, 
simplified since y of the gas сап be considered constant, and 
temperature, pressure, and density ratios are simply the 
ideal-gas functions of Mach number. 

Shock-expansion method.—This method of calculating 
supersonic flow of an ideal gas at the surface of an airfoil is 
well known, entailing simply the calculation of flow at the 
nose with the oblique-shock equations and flow downstream 
of the nose with the Prandtl-Meyer equations. Determina- 
tion of airfoil characteristics in this manner requires only a 
small amount of time, of course, compared to that involved 
when the method of characteristics is used, hence, the 
advantage of the former method. ‘The questions arise, 
however, as to exactly what the simplifying assumptions 
underlying the shock-expansion method are, and what form 
the method takes (for calculative purposes) when the gas 
displays varying vibrational heat capacities. 

The matter of simplifying assumptions may perhaps best 
be considered by employing equations. (12) and (13), the 
basic compatibility equations. If these expressions are 
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resolved into the streamline direction and combined, noting 
that 


en 2 s B др al (21) 
and 
ое Tease (30; 80 2 
there is then obtained the relation 
КЕС 
op _ 05/0C; 2үр 06 (23) 
де 1 493/00, sin 28 Os 
95/00; 
defining the gradient of р along s. If flow along streamlines 


downstream of the nose is of the simple алынады type, 
however, we have 


Op | 2vyp 06 
Os sin 28 08 (24) 


Hence, it is evident that the requirement for this type flow is 


05/00; 
95/90; 


> 1 00 (25) 


Equation (25) 1з, of course, simply an approximate statement 
of a well-known property of Prandtl-Meyer flows; namely, 
that flow inclination angles are essentially constant along 
first-family Mach lines. It follows from equation (12) that 
if equation (25) holds, then the pressures will also be еввор- 
tially constant along these lines. It does not follow, how- 
ever, that the, Mach number will be constant or, for that 
matter, that the first-family characteristic lines will be 
straight (as is the case for isentropic expansion flows about 
a corner). In fact, it may easily be shown that the Mach 
number gradient along О, is proportional to the local entropy 
gradient normal to the streamlines and that the C, lines are 
curved according to the change in M. Thus we seo that 
there is really only one basic assumption underlying the 
shock-expansion method; namely, disturbances incident on 
the nose shock (or, for that. matter, any other shock) are 
consumed almost entirely in changing the direction of the 
shock. Іп this regard, it is interesting to note that the 
assumption of Thomas (ref. 13), that pressure is a function 
only of flow deflection angle and entropy, is equivalent to 
this assumption. It follows, of course, that the most 
general solution obtainable with Thomas’ series ropresenta- 
tion of the pressure is that given by the shock-expansion 
method. 

With the assumption that all disturbances incident on tho 
shock wave are consumed, it is evident that the shock- 
expansion method provides a relatively simple means for 
calculating the whole flow field about ап airfoil, including 
the effects of shock-wave curvature. The details of such 
calculations are presented in Appendix B. In general, of 
course, the validity of this assumption and the accuracy of 
the shock-expansion method can only be checked by com- 


A STUDY OF INVISCID FLOW ABOUT AIRFOILS AT HIGH SUPERSONIC SPEEDS 343 


parison of calculations using this method with those using 
the method of characteristics. 

The shock-expansion method for a calorically imperfect, 
diatomic gas is readily deduced from the equations previ- 
ously obtained. For example, flow conditions at the leading 
edge of an airfoil can be evaluated with the.oblique-shock- 
Wave expressions (eqs. (18) through (20)) and the expression 
for conservation of energy (eq. (15)). The variation of flow 
inclination angle with pressure along the surface is then 
obtained by graphically шешш equation (24); namely, 


Ps gin 128 1 
- қы” I (26) 


where the variables y, p, and 8 are evaluated using equa- 
tions (14) through (17), employing the static temperature 
as a parameter. When extreme accuracy is not essential, 
this rather tedious calculation can be avoided, and a rela- 
tively simple algebraic solution of the flow downstream of 
the nose can be employed.* The details of this solution 
are presented in Appendix C. In the special case of flow 
at high supersonic speeds about slender airfoils, the whole 
calculation becomes particularly simple and warrants special 
attention. 

If it is assumed that the local surface slopes are small 
compared to 1 and, in addition, that the free-stream Mach 
number is large compared to 1, it follows that o and В are 
everywhere small compared to 1. In this case, equation (24) 
takes on the approximate form 


d | : 
d5|,7 7M (27) 


Furthermore, if it is assumed that y is constant at an average 
value ya for a particular flow field (this assumption appears 
reasonable since, in the temperature range up to 5,000? R, 
the change in y is less than about 15 percent), then the Mach 
number and pressure may be related by the simple expression 


Ya-1 
M=My (2) Еа (28) 
Equations (27). and (28) combine to yield the differential 


equation 
— (7a1) 


ORO 
та (2 a(2)—as | (29) 


which readily integrates (between N and S) to the form 


оо‏ “)دنجم سيم 


now denoting 


ала 
х g (Mon) | : (32) 


‘The babolated results of Noyes (ref. 14) may also prove useful in this case for Mach 
numbers up to 


= f(Moix) (31) - 


there is obtained from the oblique-shock equations, simplified 
to conform with this analysis,. 


Мору? — 1 


fsx) = — SSS (83) 
(Méca 2. г) н Mitox'—1) 
and 
__2ү„ М on? — (ya— 1) | 
Mer aa TI | (34) 
where 


Мусу هك‎ Mnt 1+ ert Мә») (35) 
With equations (30) through (35), the pressures on the sur- 
face of an airfoil may easily be obtained. In terms of pres- 


sure coefficient we have 
нв) eo 
Po N 

or 


0, M2 و‎ (Msby) Е — fM) (1 (| сей 1} (37) 


‘The advantage of these slender-airfoil expressions lies, of 
course, in their relative simplicity and, thus, the ease of 
calculation which is inherent to them. It may be noted in 
this regard that the functions f(M,ôx) and g(Msy) can be 
calculated once and for all with equations (33), (34), and (35), 
provided the variation of y, with Му is known. This 
calculation has been carried out for a constant value of y 
equal to 1.4 and for average values of y, assuming Ту--5009 
R. The results are presented in table I. 

It should also be noted that the slender-airfoil expressions 
of the shock-expansion method satisfy the hypersonic simi- 
larity law for airfoils first deduced by Tsien (ref. 15)? A 
necessary condition for the validity of these expressions is 
thus satisfied; however, the accuracy of the shock-expansion 
method, whether for slender airfoils or otherwise, remains to 
be investigated. 

METHODS FOR CALCULATING THE FLOW IN THE REGION OF THE 
„, LEADING EDGE 

As was pointed out previously, pressure disturbances 
emanating from the surface of a curved airfoil interact with, 
and thereby curve, the leading-edge shock wave. The 
geometry of this phenomenon near the leading edge of a 
convex airfoil is illustrated in figure 1. The pressure dis- 
turbances from the airfoil (expansion waves for a convex 
airfoil) travel along first-family Mach lines Q4. In addition 
to changing the inclination of the shock wave, the interaction 
between these disturbances and the shock produces another 


System of disturbances which travel along second-family 


Mach lines C; from the shock wave to the surface. 
Method of characteristics.—Àn exact solution for the sur- 
face pressure gradient and shock-wave curvature at the leading 


5 For a given value of То, Ту, to the accuracy of this analysis, is the {deal-gas function 
of Мух. Thus, knowing Ty, yx can be determined. The average value of y used 18 
a LOD Ga тд. 

6 This fact was employed by Linnell (ref. 16) to obtain an expression for pressure coefficient 
equivalent to equation (37) for the caso of constant у and to obtain explicit solutions for the 
lift, drag, and pitching-moment coefficients of several airfoils at hypersonic speeds. 


344 


Shock wave 
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FIGURE 1.—Schematic diagram of supersonic flow past a curved 
sharp-nose airfoil. 


edge may be determined in the following manner. Itis clear, 
referring to figure 1, that the change in flow angle between 
points A and C given by the compatibility equations (see eqs. 
(12) and (13)) along the path ABC must equal the change 
determined by the airfoil surface from A to С.. Similarly, 
the difference in pressure between points B and D given by 
the compatibility equations along the path BCD must equal 


that determined by the change in shock-wave inclination . 


between B and D. In reference 9, these conditions were 
employed at the leading edge to obtain equations, in simple 
parametric form, for determining the surface pressure 
gradients and shock-wave curvatures. These equations can 
be written in the form 


1 08/0C: 
1 аР | 0, | 2yP (38) 
кат” 128/00, sin 28 
05/0C 
for the surface pressure gradient, and 
| КАТ 
sin (В--о--6)-Р sin (8--о-8). 
K xi 08/00 (39) 


4 28/00; 08 
tasa. 55/9 Z) sin В 


for the shock-wave curvature, where 


COAT 


Fo) 98 


It should be pointed out that equation (38) is, of course, 
equivalent to equation (23). 


96/0C, . 
d/d | 


(40) 
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A procedure for evaluating equations (38) and (39) for a 
calorically imperfect, diatomic gas, as well as for an ideal 
gas, 18 presented іп Appendix О. Since they are exact for 
two-dimensional, steady, inviscid flow, values of the surfaco 
pressure gradient and shock-wave curvature at the leading 
edge, determined using these equations, may be put to two 
uses. First, the accuracy of approximate methods of cal- 
culating the flow field about a curved airfoil can be evaluated 
at the leading edge by comparing the values of the surface 
pressure gradient and shock-wave curvature predicted by 
these approximate methods to the values obtained using 
equations (38) and (39). Second, the pressure gradient and 
shock-wave curvature can be used to calculate the initial 
points of a characteristic solution for the flow about an airfoil. 

Shock-expansion method.—One approximate solution for 
the surface pressure gradient and shock-wave curvature has 
already been indicated in the previous discussion of the shock- 
expansion method. ‘The requirement for the application of 
the shock-expansion method is given by equation (25). 
Hence, it is apparent that under this condition the expres- 
sion for the surface pressure gradient (eq. (38)) reduces to 


Similarly, the expression for the shock-wave curvature (eq. 
(39)) reduces to 


Ж. =з 8010) o-+d) 


Ko (B) sng ‚эш B 


It should be realized that the flow field is determined by 
the basic flow equations in conjunction with the shock wave 
and airfoil surface as boundary conditions. Thus, the addi- 
tional requirement for this shock-expansion method of zero 
pressure gradient along first-family Mach lines means that 
one of the flow relations cannot be satisfied exactly (i. o., the 
flow field is overdetermined). Equations (41) and (42) | 
satisfy the shock relations and the airfoil surface as boun- 
dary conditions; however, the compatibility equations are 
only approximately satisfied. ' (See Appendix В.) 

The error in surface pressure gradient associated with 
neglecting the reflected disturbances might be expected to 
be largest in the region of the leading edge of a curved airfoil 
due to the close proximity of the shock wave and the surface. 
The magnitude of the error in this region may be deduced, 
of course, from the ratios of values of surface pressure 
gradient and shock-wave curvature given by the character- 
istics method to those given by the shock-expansion method. 
The surface~pressure-gradient ratio and shock-wave-curva- 
ture ratio can be written (using eqs. (38) and (41)) 
1 _ 08/00, 

‚ у—-—_990@ оё/дО» 
1 4 05/00; 
95/00; 
and (using eqs. (39) and (42)) 


sin (G@—o+8)+( S204 


(1 ШЕН, sin (8—0 +8) 


(42) 


(43) 


sin (8--о--6) 


(44) 


K= 
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respectively. A procedure for evaluating equations (43) and 
(44) for the flow of а calorically imperfect gas, as well as for 
an ideal gas, is presented in Appendix D. (The application 
of eq. (43) for an ideal gas has already been given in ref. 9.) 

As was stated previously, the pressure gradient and shock- 
wave curvature can be used to calculate the initial points of 
a characteristic solution for the flow about an airfoil. It is 
apparent that the initial points for a shock-expansion solu- 
tion can be found in а similar manner. With either type of 
solution, flow conditions at an initial point on the surface 
downstream of the leading edge can be calculated with the 
Ма of the appropriate value of the surface pressure gradient. 
Similarly, the flow conditions at an initial point on the shock 
wave can be obtained with the aid of the corresponding 
value of the shock-wave curvature. Additional points can 
be obtained between these two by linear interpolation. If 
the two initial points are chosen as the ends of a first-family 


Mach line, there is sufficient information available to deter- . 


mine the curvature of this Mach line. Therefore, if detailed 
knowledge of the flow in the region of the leading edge is re- 
quired, the surface, shock wave, and first-family Mach line 
can be approximated by circular arcs. (See ref. 17.) 


INVESTIGATION OF THE FLOW ABOUT AIRFOILS AND 
DISCUSSION OF RESULTS 


This study is divided into two sections: first, an investiga- 
tion of the flow in the region of the leading edge of curved 
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airfoils; and second, an investigation of the complete flow 
field about an example airfoil. Each of these sections is 
further subdivided into а consideration of the effects of 
Mach number, assuming air behaves as an ideal gas, and 
into a consideration of the combined effects of Mach number 
and gaseous imperfections. In the latter regard, the prin- 
cipal emphasis is placed on the caloric imperfections pre- 
viously discussed. 


FLOW IN THE REGION OF THE LEADING EDGE OF CURVED AIRFOILS 


Ideal-gas flow.—The results of the calculations (using 
eq. (D4)) of the surface pressure gradient are presented in 
table П and figure 2.” The values presented in the table are 
for a range of Mach numbers from 1.5 to о and for leading- 
edge deflection angles from 0° to 45°. Where no value 
appears in the table, the flow behind the shock wave is sub- 
sonic. ` Corresponding results of the calculations of surface- 
pressure-gradient ratio are presented in table ПІ and figure 3. 
From these results it is seen that except near shock detach- 
ment, the surface-pressure-gradient ratio varies only from 
0.98 to 1.02 for Mach numbers less than 4. Therefore, very 
little error will result from the use of the shock-expansion 
method for the surface pressure gradient at these lower 
Mach numbers. For Mach numbers greater than 4 (even 

1 Oharts were also presented for surface pressure gradient, surface-pressure-gradient ratio, 
and shock-wave curvature for ideal-gas flows in reference 9; however, the results given in the 


present report are somewhat more extensive. These results were also presented in cross- 
plotted forms in reference 17. 
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FIGURE 3.— Variation of surface-pressure-gradient ratio with deflection 
angle for various free-stream Mach numbers. 


including о), the ratio only varies between about 0.9 and 
1.1 (again except near detachment) and, therefore, only small 
error might be expected from the use of the shock-expansion 
method. Near detachment, at the higher Mach numbers, 
the surface-pressure-gradient ratio attains a very large range 
of values and the maximum value increases with Mach num- 
ber. For these conditions, then, the use of the shock-ex- 
pansion method for calculating the surface pressure gradient 
at the leading edge would result in appreciable error. 

The flow along the surface is isentropic. Hence, it can be 
shown that y, the ratio of surface pressure gradient to that 
given by the generalized shock-expansion method, is also the 
velocity-gradient ratio, the Mach number gradient ratio, the 
Mach angle gradient ratio, the density-gradient ratio, and 
the temperature-gradient ratio. Any of these gradients may 
be found, then, by calculating the gradient, using the shock- 
expension method and applying the appropriate value of y. 
This property of the ratio y makes it useful in the application 
of the method of characteristics with any of the coordinate 
systems commonly employed in the compatibility equations. 

The results of the calculations of shock-wave curvature 


% For the particular case of infinite free-stream Mach number and zero deflection angle, the 
pressure-gradient ratio 1s double valued. From equation (D8), it is apparent that y ts unity 


3 / T(v—1) 
zi V 3~ as the 


ог zero deflection. Yet, at infinite Mach number. ў approaches 25 
deficction angle approaches тего. 
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(eq. (D5)) are presented in table IV and figure 4.° Similarly, 
the results of the calculations of shock-wave-curvature ratio 
(eq. (44)) are presented in table V and figure 5. Except near 
detachment, the curvature ratio varies from 0.92 to 1.08 for 
all Mach numbers including o. Thus, only small errors 
would result from using the value of the shock-wave curva- 
ture given by the shock-expansion method for all flow con- 
ditions except near shock detachment. 
Calorically-imperfect-gas flows.—With increasing Mach 
number and leading-edge slope, the temperature ratio across 
an oblique shock wave increases as shown in figure 6. Аз 


· the temperature behind the shock wave increases, the be- 


havior of air diverges from that of an.ideal gas as discussed 
previously. Below 800° R., the divergence is not significant, 
and the equations for ideal-gas flow can be applied with only 
minute errors resulting. Above 800° R., the energy of the 
vibrational degrees of freedom of the gas molecules is appre- 


-ciable and becomes greater with increasing temperature. 


For these conditions, the specific heats and their ratio vary, 
significantly with temperature. The equations developed 
previously consider these effects and permit the extension of 
the solution for surface pressure gradient and shock-wave 
curvature to the case of calorically imperfect gases. These 


3.0 HHH 
1 a B 


3.5 


Ke 
Ku 


Shock-wave curvature 


Leading-edge deflection angle, 8, degrees 


Ғісонв 4.—Variation of shock-wave curvature with leading-edge 
deflection angle for various free-stream Mach numbers. 
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FIGURE 6.— Variation of leading-edge deflection angle with free-stream 
Mach number for various values of the temperature ratio. 
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equations are valid for temperatures up to the order of 5,000° 
R. For a free-stream temperature of 500° R., therefore, the 
shaded area between lines of constant temperature ratio of 
1.5 and 10, in figure 6, represents the range of conditions for 
which the method developed in this report for the flow of a 
diatomic, calorically imperfect gas would apply. 

The excitation of the vibrational degrees of freedom of the 
gas molecules requires а finite number of collisions, causing 
the well-known heat-capacity lag discussed in references 
5 and 6. The flow distance (1. e., along the streamline) 
required to establish equilibrium is usually small in dense air 
and will be considered infinitesimal in this report. Also, the 
dissociation of air (see ref. 6) will not be considered here. 

Since the free-stream static temperature is an additional 
parameter in calculations of flow of imperfect gases, only a 
limited number of calculations of (1/K..) (GP/dW), y, К./К,, 
and к were made. Тһе purpose of these calculations is to 
compare the variations of these quantities with the values 
as given by the ideal-gas-flow computations. The calcula- 
tions followed the procedure described in Appendix D. A 
free-stream static temperature of 500° R. was used. The 
results of these calculations are presented in table VI for 
various Mach numbers and leading-edge deflection angles. 

The surface pressure gradients for an ideal gas and for а 


‘calorically imperfect, diatomic gas are compared in figure 7. 


In all cases calculated, the gradient is smaller for the im- 
perfect gas and diverges gradually, with increasing free- 
stream Mach number and deflection angle, from the value 
of the gradient for an ideal gas. This divergence is consistent 
with the increasing effects of the caloric imperfections due 
to the increasing temperature behind the shock wave. 

The surface-pressure-gradient ratio for the imperfect gas 
is compared in figure 8 with the ratio for an ideal gas. The 
surface-pressure-gradient ratio is smaller for impurfect-gas 
flows than for ideal-gas flows indicating that the effects of 
shock-wave and expansion-wave interaction are greater. 
This result is attributed, in part, to the fact that the angle 
between the shock wave and the airfoil surface is smaller for 
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FIGURE 7.—Comparison ‘of the ‘variation of surface pressure gradient 
with leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 
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numbers for an ideal gas and for a calorically imperfect, diatomic gas. 
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ratio with leading-edge deflection angle for various free-stream Mach 
numbers for an ideal gas and for a calorically imperfect, diatomic gas. 


differ only by about 5 percent (see fig. 8). 


_to disturbances incident on the wave. 
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imperfect-gas flows than for ideal-gas flows. The difference 
between the imperfect- and ideal-gas calculations increases 
with increasing Mach number and deflection angle, as might 
be expected. In the range of Mach numbers and flow deflec- 
tion angles investigated, however, the extreme values of y 
Thus, it is appar- 
ent that while the shock-expansion method will not be quite 
as accurate for calorically-imperfect-gas flows as for ideal-gas 
flows, the method will not be expected to be invalidated. 

A divergence with Mach number and deflection angle is 
also apparent in figure 9 in which the shock-wave curvatures 
for an ideal gas and a calorically imperfect, diatomic gas are 
compared. This divergence is compatible with the change 
in surface pressure gradient due to the caloric imperfections 
of the gas. The shock-wave-curvature ratio for a calorically 
imperfect, diatomic gas and this ratio for an ideal gas are 
shown in figure 10. Again, it is seen that the effect of caloric 
imperfection is to increase the effects of shock-wave and 


expansion-wave interaction. 


COMPLETE FLOW FIELDS 


Ideal-gas flows.—The effects of Mach number of primary 
interest here are, of course, those which result from the inter- 
action between the leading-edge (or other) shock wave and 
small disturbances originating on the surface of an airfoil. 
Although the reflected disturbances that are the product of 
this interaction will have the largest effect on the flow near 
the leading edge, their influence on the complete flow field 
about an airfoil also warrants investigation. Further in- 
sight into these effects can be obtained in the region just 
downstream of the shock wave without regard for the shape 
of the airfoil producing the shock. То this end, it is con- 
venient to consider the ratio Raa NA орде (see eq. (25)) 
which may be termed “the disturbance strength ratio” since, 
in the region under consideration, it is а measure of the ratio 
of strengths of disturbances reflected from the shock wave 
This ratio may be 
evaluated with equation (40). This calculation has been 
carried out for Mach numbers from 3.5 to о (y=1.4) and 
flow deflection angles approaching those corresponding to 
shock detachment (1. e., М, == 1), and the results are pre- 
sented in figure 11. It is evident that except near М,=1, 
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the ratio is small (in absolute value) compared to 1 through- 
out the entire range considered—this observation also ap- 
' plies, of course, at lower supersonic Mach numbers. Thus 
it 1s indicated that almost all of an incident disturbance is 
generally absorbed in the shock wave, provided the air be- 
haves like an ideal diatomic gas. This result is substan- 
tially the same, of course, as that which is assumed in deriv- 
ing the shock-expansion method of calculating flows about 
airfoils, and therefore yields additional credencein the method 
for high Mach number as well as low Mach number applica- 
tions. lt should also be noted that this result is contrary 
to that obtained by Lighthill (ref. 18), who reports that for 
hypersonic Йозув(А% 6>>1) a disturbance is reflected from 
а shock wave with opposite sign but essentially undiminished 
strength, Lighthill's conclusion is based on an incorrect 
evaluation of his results for the case of very high Mach 
numbers. 

As an over-all check on the shock-expansion Method: 
surface pressure distributions calculated thereby are com- 
pared in figure 12 with those obtained with the method of 
characteristics for a 10-percent-thick parabolic-arc biconvex 
airfoil (a—0?) operating at free-stream Mach numbers of 
3.5, 10, and о. (Additional calculations presented in ref. 
19 were also performed for Мо=5, 7.5, and 15.) Predictions 
of the slender-airfoil approximation to the former method 
for high supersonic speeds are also shown. There is no 
apparent difference between the pressure distributions given 
by the method of characteristics and the shock-expansion 
method at a Mach number of 3.5; at a Mach number of 10, 
nnd more so at infinite Mach number, however, the latter 
method predicts pressures which are slightly low downstream 
of the nose. This result would be deduced from figure 11 
where it is observed that, at the Mach numbers under con- 
sideration, expansion waves incident on the nose shock wave 
are reflected back toward the surface as compression waves 
of relatively small but increasing strength with increasing 
Mach number. The effect of these waves does not become 
pronounced even at infinite Mach number (see fig. 12 (c)) 
and it can be seen, upon comparison of these results with 
those presented previously for the pressure-gradient ratio 
(see fig. 3) that the effect of the reflected waves dissipates 
somewhat downstream of the nose. The shock-expansion 
method is thus further substantiated as being a reliable 
simplified method for predicting the flow about airfoils at 
high supersonic speeds, again, so long as the air behayes as 
ал ideal diatomic gas. The further simplified slender-airfoil 
method also appears to be a good approximation over the 
entire range of Mach numbers," although, as would be ex- 
pected from the assumptions made in its development, it is 
in somewhat greater error than the shock-expansion method 
at lower Mach’ numbers. 

The shape of the shock wave given by the shock-expansion 
method, as presented in Appendix B, is compared in figure 
13 with the shape given by the method of characteristics 
for the biconvex airfoil at 1М-- о. (These shock waves 
correspond to the pressure distributions given in fig. 12 (c).) 
The shock-expansion method gives a reasonably good 

и The hybrid expresston for pressure coefficlent obtained by Ivey and Oline (ref. 5) gives 


reasonably good results also, although not as accurate as the slender-airfoll method at the 
higher Mach numbers under consideration. 
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approximation to the shock-wave shape though, as would 
be expected from the results given in figure 5, the curvature 
is somewhat too small. (A procedure for getting a closer 
approximation to the shock-wave shape is also given in 
Appendix B.) It can be seen, however, that this method of 
determining the shock-wave shape is far better than the 
assumption of a straight shock wave that is often associated 
with the shock-expansion method. Evidently, then, the 
shock-expansion method can also be used to calculate the 
flow in regions away from the airfoil surface. (See Appendix 
B.) 
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airfoil section at «—0? and Mi= о. 


The relative accuracy at high Mach numbers of the 
slender-airfoil method and linear and second-order potential 
theories may be seen in figure 14. As might be expected, the 
slender-airfoil method is more accurate than linear theory 
at both 14,=5 and 15 and more accurate than second-order 
theory at 140-15. Itis perhaps surprising to note, however, 
that at the lower Mach number of 5, the slender-airfoil 
method is also somewhat superior to the second-order theory. 

The pressure distributions of figures 12 and 14 have been 
employed to calculate the zero-lift drag of the biconvex 
airfoil, and the results of these calculations, along with 
additional predictions of the different methods, are shown 
in figure 15. Predictions of the shock-expansion method are, 

_of course, in best agreement with those of the method of 
characteristics; while the slender-airfoil method, although 
slightly less accurate than the shock-expansion method, is 
apparently superior to both linear and second-order theories 
at Mach numbers above 3.5. 

The preceding findings verify that, so long as the dis- 
turbance strength ratio is small compared to 1, the flow 
along streamlines is essentially of the Prandtl-Meyer type. 
If we choose, on the basis of these findings, a maximum 
absolute value for 535 Cı of 0.06 (note the maximum value of 

6/0 Cs 


i A for the cases presented in fig. 12 was approximately 
0.06 at I= о), the region in which the shock-expansion 
method is applicable can readily be obtained from figure 11. 
The upper boundary line of this region is shown in figure 16, 
and it is evident that it lies only slightly below (about 1°, 
in general) the line corresponding to shock detachment given 
approximately by the M.=1.0 line. Almost the entire 
region of completely supersonic (ideal gas) flow is then 
covered by the method. (See shaded area of fig. 16.) The 
range of applicability shown in figure 16 is appreciably 
larger than that indicated by Rand (ref. 20) who required 


Pressure coefficient, Со 


REPORT 1123—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


с 
© 

Ax 
c 
V 
o 
“- 
че 
يو‎ 
о 
o 
Ф 
ъъ 
ә 
يو‎ 
„ә? 
о» 
| 
а 


Ш 
= 
БЕ 
^ 


NSE 


O 


Method of characteristics 
—-— Slender-airfoll method 
-------- Linear theory 
—— — Second-order theory 


190 р 6 ‚8 LO 
Chordwise station, А/С 
(а) Ме==5 
(b) Mo=15 


FIGURE 14.—Pressure distribution on 10-percent-thick biconvex airfoil 
| section at «=0°, . 


that the entire flow field be of the true Prandtl-Meyer type 
(i. e., that all flow properties be constant along first-family 
Mach lines and not just 6 and p). The results presented in 
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FIGURE 15.—Variation of drag coefficient with free-stream Mach. 
number for 10-percent-thick biconvex airfoil section at a=0° (y=1.4). 
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figures 12 and 13 show, however, that this restriction is not 
necessary. 

. The question naturally arises concerning the corresponding 
range of applicability of the slender-airfoil method. This 
question may be answered, 1п part, by comparing separately 
the predictions of the method for oblique-shock flows and 
expansion flows with those of the exact oblique-shock 


equations and Prandtl-Meyer equations. Such а com- 


parison is shown in figure 17 in terms of the percentage 
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error in the pressure coefficients predicted by the slender- 
airfoil method. As would-be expected, this method does 
not exhibit good accuracy over the wide range of appli- 
cability of the shock-expansion method; however, it is 
indicated that it should predict pressure coefficients with 
less than 10-percent error down to Mach numbers as low 
as 3 for airfoils producing flow deflections up to as high as 25°. 

As a further check on the utility of the slender-airfoil 
method, the pressure coefficients on the 10-percent-thick 
biconvex airfoil have been calculated with this method and 
the shock-expansion method at a Mach number of 10 and 
angles of attack up to about 30°.% The results of this 
calculation are shown in figure 18 (see fig. 12 (b) for a=0°) 

1? These conditions are within the range of applicability of the shock-expansion method 
as defined in figure 16; hence, the use of the method as a base of comparison seems justified. 
Віпсе the shock-expansion method 15 far less tedious to apply than the method of character- 


istics, it will be employed as such a base іп subsequent calculations whenever the conditions 
teing investigated have been determined to be within its range of applicability. 


» ? 


Ss 
= 
S 
A 
D 
o 
о 
5 
Calorically imperfect gas, 7,=500°R 
e Shock-expansion method 
= ------- Slender-oirfoil method утуд 
© Ideal gas 
С; Shock-exponsion method ууу 
D 
o 
o 
сл 
2 
о 
$ 
= 
ف‎ 
5 
Ет 
о 
o 
© 
نه‎ 
Е 
о 
= 


Angle of attack, a, degrees 


(a) Lift coefficient. 
(b) Drag coefficient. 


(c) Moment coefficient. 


FraunE 19.—Variation of force and moment coefficients with dais of 
attack for 10-percent-thick biconvex airfoil section at M,—10 


(ү--1.4). 


where it is seen that the agreement is reasonably good, even 
at the highest angle of attack. This fact is reflected іп 
figure 19 showing the force and moment coefficients for the 
airfoil as a function of angle of attack. Little difference 
is observed in the force coefficients as calculated by the 
two methods, while the moment coefficients display more 
pronounced but, nevertheless, small differences at the 
higher angles of attack. 

From these and previous considerations, the ranges of 
applicability of the shock-expansion and slender-airfoil 
methods for supersonic ideal-gas flows are reasonably well 
established. It remains now to determine the manner and 
extent to which gaseous imperfections in the flow at higher 
supersonic speeds may alter these ranges, and the reasons 
therefor. | 

Imperfect-gas flows.—As a first step toward investigating 
the effects of gaseous imperfections on the high Mach 
number flows under consideration, it 18 conyenient to extend 


our consideration of the disturbance strength ratio RE 
2 


It is recalled that when air exhibits a constant value of y. 


equal to 1.4 (the value for an ideal diatomic gas), the disturb- 
ance strength ratio is small at arbitrarily large Mach numa- 
bers, provided the flow deflection angles are not too close 
to those for shock detachment. One of the most important 
effects of gaseous imperfections is, however, to decrease y 
of the disturbed air below this value due to the excitation of 
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additional degrees of freedom (e. g., vibrational) in the 
molecules at the high temperatures encountered at high 
Mach numbers. Indeed, at arbitrarily high Mach numbers 
it might be expected that y of the disturbed air would 
approach 1, since the number of degrees of freedom p 
effectively become very large (see, e. g., refs. 3 and 6). 

this case, however, the extent of the disturbance flow = 
is decreased {о а layer at the surface of the body which is 
negligibly thin compared to that for the case of ideal-gas 
flow. Thus, it is apparent that significant changes in the 
flow about airfoils at high Mach numbers may result from 
decreases in y of the disturbed air; hence, the effects of such 
decreases on the disturbance strength ratio would appear to 
warrant attention. 

A detailed analysis of these effects is impractical at the 
present time, due to the limited range over which the vari- 
ation of y with temperature is accurately known. Even in 
the range where this variation is so known, the additional 
complication required to consider the effects of variable y 
and the addition of another independent parameter (free- 
stream temperature) make extensive calculations of the 
disturbance strength ratio impractical. However, some 
knowledge of these effects can be gained by performing the 
calculations for one free-stream Mach number and tem- 
perature. Such calculations have been carried out at o 
Mach number of 10 for a free-stream temperature of 500° 
R. and the results are presented in figure 20. Тһе curve for a 
calorically imperfect gas cannot be extended to shock 
detachment because the temperature behind the shock wave 
exceeds that for which the calorically-imperfect-gas equations 
are valid. It can be seen that the effect of the caloric im- 
perfection of air is to increase the value of the disturbance 
strength ratio and that the effect increases with increasing 
temperature or decreasing y. However, it appears that if y 
does not decrease appreciably below 1.3, as in this case, the 
disturbance strength ratio is still small compared. to unity. 
It might be expected, therefore (as previously found for flow 
in the region of the leading edge), that the shock-expansion 
method would continue to predict the flow about complete 
airfoils with reasonable accuracy. This pomt has been 
checked with the methods developed previously for analyzing 
the flow of a calorically imperfect, diatomic gas at local air 
temperatures up to about 5,0009 К. (note y has a value only 
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FIGURE 20.—Effect of the caloric imperfections of air on the disturbance 
strength ratio at М0=10 and Ty—500? R.* 
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FIGURE 21,—Pressure distribution on lower surface of 10-percent-thick 
biconvex airfoil section at My=10, 71-5009 R, а=19.9°. 


slightly less than 1.3 at this temperature). In particular, the 
pressure distribution on the lower surface of the biconvex 
airfoil at Mo=10, а--19.99 and Тұ--5009 К. (Ту عع‎ 40009 К.) 
has been calculated with both the method of characteristics 
and the shock-expansion method.? The results of these 
calculations are presented in figure 21, andit would appear that 
the conclusions drawn from figure 20 aresubstantiated. (Due 
to the lower temperatures, pressures in the expansion flow 
about the upper surface are not influenced by caloric imper- 
fections and, hence, are the same as shown in figure 18 (b).) 
For these same conditions, the shapes of the shock waves 
given by the shock-expansion method and by the method of 
characteristics (both for a calorically imperfect gas) are com- 
pared in figure 22. Just as in the case of ideal-gas flows, the 
shock-expansion method gives a good approximation to the 
shock-wave shape, far better than the assumption of a 
straight shock wave. Thus, it is seen that the conclusion 
drawn from figure 20 should also apply for the use of the shock- 
expansion method to calculate the flow field away from the 
airfoil surface. 

Shown also in figure 21 is the pressure distribution ob- 
tained by the shock-expansion method for an ideal gas 
(yi—1.4). Itis apparent, on comparing this pressure distri- 
bution with the other distributions, that although the effect 
of caloric imperfections on the disturbance strength ratio 
is small, the pressures are appreciably reduced by the increase 
in specific heats. The extent of this reduction is more com- 
pletely illustrated in figure 23 where the lower-surface pres- 
sure distributions on the biconvex airfoil are presented for 
M,=10 and Т,--5009 R., at «—=0°, 10°, 19.9°, and 30°. Ав 
one might expect, the reduction in pressures increases with 
angle of attack (due to the corresponding increase in static 
temperature of the disturbed air). The pressure coefficients 
calculated with consideration of the imperfections in the gas 
are less on the lower surface (up to 6 percent at the leading 
edge and 15 percent at the trailing edge) than those calculated 
assuming the gas behaves ideally. The upper-surface pres- 

и For added enso of саспа оп the expansion method of Appendix О vas employed. This 
method is also employed іп ali subsequent calculations of this type since it has been found to 


yield results differing by less than 1 percent from those obtained by the more tedious graphical 
integration method. 


821606--б0----24 


(up to about 10 percent for а--309). 
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sures are again unaffected by the caloric imperfections of air 
in all the cases presented (except at a=0°) since this surface 
experiences lower pressures and, hence, lower temperatures. 
They are therefore the same as shown in figure 18. Shown 
also in figure 23 are the pressure distributions calculated 
with the slender-airfoil method for y=ya. The accuracy of 
this simplified method is substantially the same as was pre- 
viously observed for the corresponding method in the case 
of ideal-gas flows, although the local error may be greater 
iban the reduction in pressure coefficients due to the caloric 
imperfections of air. 

The force and moment coefficients, corresponding to the 
lower-surface pressure distributions shown in figure 23 and 
the upper-surface distributions of figure 18, are presented in 
figure 24. The reduction in the lower-surface pressures 
leads, of course, to а general reduction in all three coefficients 
The slender-airfoil 
method again predicts these coefficients with surprising - 
accuracy. | 

In order to further assess the accuracy of the slender-airfoil 
method, some additional calculations were carried out for the 
biconvex airfoil аф a=0° and 14-20 and 30. The pressure 
distributions for these cases were calculated by the shock- 
expansion method, slender-airfol method (у=), and 
slender-airfoil method (у==ү;). These results are presented 
in figure 25, and it is observed that the use of y, rather than | 
y. improves the accuracy of the slender-airfoil method. 
The extent of this improvement in the case of drag coefficient 
is shown in figure 26; it would appear that predictions of the 
slender-airfoil method (ү--та) and shock-expansion method 
are in as good agreement as for ideal-gas flows (see fig. 15). 
On the basis of these and previous results, it may be con- 
cluded, then, that not only does the shock-expansion method 
retain its range of applicability when air exhibits caloric 
imperfections, provided y of the disturbed air is not appre- 
ci&bly less than 1.3, but also the slender-airfoil method 
(у= үл) retains its range of applicability. 

It would be surprising indeed, however, if this conclusion . 
continued to apply as y of the disturbed fluid approached 1 
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FiguRE 22.—Shape of the shock wave for 10-percent-thick biconvex 
airfoil section at «=19.9°, 240-10, 7T, 500? R. 
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° FIGURE 23.— Pressure distribution on lower surface of 10-percent-thick biconvex airfoil section for various angles of attack at Мо= 10 and To= 500° R. 


25/00, . 


since, as indicated previously, 35/0 increases with decreas- 


їр у. Although the manner in which у varies with tempera- 
ture is not known in this range, some knowledge of these 
effects can be gained by repeating the ideal-gas calculations 
for constant values of y between 1.4 and 1.0.4 Such calcüla- 
tions have been carried out for infinite Mach number since, 
05/01 
95/06, 
and the results are presented in figure 27. It is seen that 
except near detachment, the disturbance strength ratio ш- 
creases with decreasing y. ‘This increase is slow at first; 
spre is still less than 0.1 at 
y=1.3. This result is in agreement with the previous con- 


in this case, has its maximum value for a given y, 


for example, the value of 


и Since the enthalpy is negligibly small compared to the mass kinetic energy of the undis- 
turbed fluid at the high Mach numbers of interest and, hence, у of this fluid does not influence 
the flow, this approach corresponds to employing an average value of у for the disturbed 
fluid, 


clusion regarding cases where y is greater than 1.3. How- 
02/20: 
05/0C, 
fact, approaches 1 аз у аррговсһев 1. Тһе effect on pressure 
distributions of this increase in the strength of the reflected 
disturbances may also be investigated by using the ideal-gas 
relationships in combination with appropriate values of y. 
The limiting case of infinite free-stream Mach number and 
y=1.0 (for the disturbed fluid, see footnote 14) has already 
been investigated by Busemann (ref. 21) and more recently 
by Ivey, Klunker, and Bowen (ref. 22). Іп this case, as 
pointed out previously, the shock wave emanating from the 
leading edge remains attached to the surface downstream of 
the leading edge, (this is easily verified with the oblique- 
shock-wave equations), and the disturbance flow field is 
confined to an infinitesimally thin layer adjacent to the 
surface. In addition, the velocity along a streamline down- 
stream of the shock is constant, as may easily be shown with 


ever, continues to increase as y decreases, and, in 
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FIGURE 24.—Variation of force and moment coefficients with angle of 
attack for 10-percent-thick biconvex airfoil section at |Мо--10, 
Т 07-5009 R, 


the compatibility equations. Surface pressures therefore 
become a simple function of airfoil geometry 


Q,—2 sin? 081-2 COS Óg os | sin ба dx 
0 


5 (45) 


varying, to a firs& approximation, directly with the square 
of the component of free-stream velocity normal to the sur- 
face (i. e., the flow is approximately of the Newtonian 
corpuscular type). With this theory, then, and the method 
of characteristics, we can,get an idea of both the extent to 
which changes of y from 1.4 toward 1 will alter surface 


pressures, and the accuracy with which the shock-expansion ' 


theory predicts the alterations. To this end, figure 28 is 
presented showing the pressure distributions about the 
biconvex airfoil at Му= о as calculated by the several 
methods for different values of y. It is observed that, 
Whereas the shock-expansion method agrees very closely 
with the method of characteristics for y=1.4, there is a 
large difference at y=1.05. This, of course, is precisely 
what one would expect from the previous discussion of the 
disturbance strength ratio. On the other hand, if the two 
characteristics solutions and the Busemann method are 
considered in order of decreasing y, it is indicated that the 
characteristics solutions approach the Busemann theory as 


ж 
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ү approaches 1. Fory=1.0 and M= œ the shock-expansion 
method, in turn, predicts & discontinuous pressure distribu- 
tion with a pressure coefficient equal to that of the Busemann 
theory at the leading edge but & pressure coefficient of zero 
at all points downstream of the leading edge. Hence, it 
may be concluded that when the free-stream Mach number 
approaches infinity and у approaches 1, the Busemann 
method rather than the shock-expansion method for calcu- 
lating the flow about airfoils should be employed. 
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FIGURE 25.—Pressure distribution on 10-percent-thick biconvex airfoil 
section аб «--09 and 7525500? R. 
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FIGURE 26.—Variation of drag coefficient with Mach number for 
10-percent-thick biconvex airfoil section for «—0? and 75:500? R. 
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FIGURE 27.—Variation of disturbance strength ratio with deflection 
angle at infinite free-stream Mach number for various values of w. 


CONCLUSIONS 


Inviscid flow about curved airfoils at high supersonic 
speeds was investigated analytically, first assuming air be- 
haves as an ideal gas, and then assuming it behaves as a 
thermally perfect, calorically imperfect gas. This study has 
led to the following conclusions: 

1. So long as air behaves as an ideal gas, the shock- 
expansion method may be used with good accuracy to pre- 
dict the flow about a curved airfoil up to arbitrarily high 
Mach numbers, provided the flow deflection angles are 
about 1° or more below those corresponding to shock detach- 
ment. ‘This conclusion applies not only to the determina- 
tion of surface pressure distributions, but also to the deter- 
mination of the whole flow field about the airfoil. 

2. An approximation to the shock-expansion method, 
applicable to ideal-gas flows about slender airfoils at high 
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Еісонв 28.—Pressure distribution on 10-percent-thick biconvex airfoll 
section at М= o and a=0°. 


Mach numbers, predicts pressure coefficients in error by less 
than 10 percent for Mach numbers above 3 and flow deflec- 


tion angles up to 25°. ы 

3. So long as caloric imperfections of air do not decrease 
the ratio of specific heats appreciably below 1.3 (correspond- ' 
ing to air temperatures up to the order of 5,000? R.), the 
shock-expansion method, generalized to include the effects 
of these imperfections, is substantially as accurate as for 
ideal-gas flows. The principal effect of caloric imperfections 
is to reduce pressure coefficients by as much as 15 percent. 
The slender-airfoil method can also be made as accurate as 
for ideal-gas flows by employing an average value of the 
ratio of specific heats. 

4. If the ratio of specific heats approaches 1, as it may at 
extremely high Mach numbers, the shock-expansion method 
саһ be in considerable error. In this case, the Busemann 
method for flow in the limit of infinite Mach number and 
specific-heat ratio of 1 applies with reasonable accuracy. 


АмЕв AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY Coamurren FOR ÁERONAUTICS, 
Morrstr FELD, Caur., January 9, 1952. 


APPENDIX A 


METHOD OF CHARACTERISTICS FOR TWO-DIMENSIONAL 
FLOW OF A CALORICALLY IMPERFECT GAS 


In the application of the method of characteristics for a 
calorically imperfect, diatomic gas to the particular problem 
of analyzing the flow about curved two-dimensional airfoils, 
many of the calculations are identical to those encountered 
in the solution of any problem where characteristics theory 
is employed. Since the details of these calculations are well 
known. and well reported (see, e. g., ref. 12), they will not be 
repeated here. 

A lattce-point system with an initial-value, numerical 
computing procedure will be used. The form of the com- 
patibility equations to be employed was developed pre- 
viously; however, it is convenient for purposes of calcula- 
tion to substitute the pressure ratio, 2/4 into these equations 
ond to rewrite them as difference equations. Equations 
(12) and (13) are thus reduced to the following forms: 


(р/ф)в— (Plo) a= — da (50— 84) 
and 
(т/аде-- )2/00( s— a(8e— ôn) (А2) 
where ко 
yY pIa 
ACTAS. (A3) 


lt is also convenient to employ several reference curves. 
These curves can be divided into two groups. The general 
reference curves consist of y and A(T) as a function of 
temperature, 7. Equations (14).and (17) are used to deter- 
mine these curves. A second set of shock-wave reference 
curves consisting of p/q, 6, and 6 as a function of tempera- 
ture, 7, are determined by use of equations (18) through 
(20)—the values of Т ала My, are presumed known. 

In the computations three types of points are encountered. 
These are (1) & point in the flow field between the shock 
wave and the airfoil surface, (2) a point on the airfoil sur- 
face, and (3) a point just downstream of the shock wave. 
Each one of these types of points requires a slightly different 
computing procedure and they will be considered in order. 


POINT IN THE FLOW FIELD BETWEEN THE SHOCK WAVE AND THE 
AIRFOIL SURFACE 


Figure 29 (a) shows a schematic diagram of the system of 
points to be considered in these calculations. Point C is 
the unknown point at the intersection of the first-family 
characteristic line passing through point A and the second- 
family characteristic line passing through ‘point B. Six 
quantities are known at both points A and B, and the prob- 
lem is to calculate these same quantities at point C. These 


1$ This form of the compatibility equations (іп p and 4 coordinates) was also used іп ob- 
taining some of the characteristics solutions for {deal-gas flows. The majority of these solu- 
tions were carried out, however, with the compatibility equations in 8,4, and entropy coordi- 
nates, since 16 was found that greater accuracy was usually obtained for в given net size. 
In general, the net size employed yielded pressures at from 30 to 35 surface points оп an airfoil 
with a maximum error In tho corresponding pressure coefficients equal to less than 1 percent of 
the pressure coefficient at tho leading edge. 


(A1) 


quantities are т, y, à, p/q, T, and T,. The first five quantities 
are of obvious significance. ‘The sixth, 7», is defined as the 
static temperature, just downstream of the shock wave, on 
the streamline passing through the point C. 

The physical coordinates of the point C(zc, yc) may be 
determined by standard procedures such as those given in 
reference 12. Іп order to determine the quantity бе, it is 
necessary to solve equations (А1) and (А2) simultaneously, 
thus, | 

js ХлбА-ҒАвбь--(р/ада--(р/4дв 
ha ths 


Equation (А1) or (А2) is then used to obtain (p/qo)c. 
There remains only the problem of determining То and 
Ts, at point C. The temperature T., is obviously constant 


(A4) 
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(a) Point in field. 
(b) Point on surface. 
(c) Point on shock wave. 
FicunE 29.— Diagram of point system in the method of characteristics 
. * for the two-dimensional flow of a calorically imperfect gas. 
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along the streamline through С. This quantity may there- 
fore be calculated in the same manner as the entropy is 
calculated in similar flow fields for ideal-gas processes (see, 
e. g., ref. 12). Furthermore, since the flow along stream- 
lines downstream of the shock wave is isentropic, equation 
(16) may be applied in the form 


(р/аде _ АЧ) 
(0/9). АС) 


The pressure, (2/9) г, is defined in a manner analogous to 
Т. and may thus be determined using the shock-wave 
reference curves and the known value of Т... Similarly, 
A(T.,) may be determined from the general reference curves. 
The only unknown in equation (А5) then is А(70) which 
may now be calculated. Once A(To) is determined, Тс may 
be determined by again using the general reference curves. 
All six quantities, zc, Ус, бс, (p/qo)e, Tc, and 7, have now 
been determined. 


POINT ON.THE AIRFOIL SURFACE 


(А5) 


Figure 29(b) shows a schematic diagram of the points to 
be considered in these calculations. "The physical coordi- 
nates of point C(zc, yc) are first calculated by solving simul- 
taneously the equation of the second-family Mach line 
passing through point B and the equation of the airfoil 
surface. When хс and yc have been determined, ёс is readily 
obtained from the equation of the airfoil surface. Equation 
(A2) is then applied to determine (p/qoc. 

Since the airfoil surface is a streamline, T5, is constant 
along the surface and may be evaluated at the leading edge. 
The temperature, Te, may then be determined using equation 
(А5) and the previously described procedure. АП six quanti- 
ties, тє, Ус, (p/qo)c, ёо, Іс, and T.,, are thus determined. 

In the special case of the first point on the airfoil surface 
downstream of the leading edge, the pressure ratio is calcu- 
lated using the pressure gradient evaluated at the leading 
edge. (See section Methods for Calculating the Flow in 


* 
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the Region of the Leading Edge, and Appendix D.) Addi- 
tional initiel points in this region may be calculated by the 
procedure previously described. 


POINT ON THE SHOCK WAVE 


Figure 29 (c) shows & schematic diagram of the points to 
be considered in these calculations. The physical coordi- 
nates of point C(xc, yc) are first calculated by solving simul- 
taneously the equation of the first-family Mach line passing 
through point A and the equation of the shock wave linear- 
ized at point D, the last known point on the wave. The 
variation of 2/00 with à along the shock wave may be approxi- 
mated by the relation 


@lade—(lado=| HE" | фе—%ь) ` (A9 


In this equation Б la | is the rate of change of р/т 


with 5 along the downstream side of the shock wave evalu- 
ated at point D. Because of the complicated nature of the 
shock-wave equations, it is generally efsier to evaluate 


E graphically or numerically from the shock-wave 


reference curves; however, this derivative may also be eval- 
uated from the equations given in Appendix D. Equations 
(А1) and (Аб) are solved simultaneously for бе, thus, 


12929] 


when ёс has been calculated, Тс and, in turn, (p/qo)o may be 
determined from the shock-wave reference curves. Since 
point C in this case is just downstream of the shock wave, 
То ала To, are identical. The.six quantities, zc, yo, (p/qo)o, 
óc, Tc, and Teg have now been determined. 


бо== 


APPENDIX B 


SHOCK-EXPANSION METHOD FOR CALCULATING THE FLOW ` 


FIELD ABOUT AN AIRFOIL 


An initial-value procedure which is similar to, although 
markedly simpler than, that associated with the method of 
characteristics may be employed to carry out this calcula- 
tion. To illustrate, consider figure 30. With the oblique- 
shock-wave and expansion equations, all fluid properties 
at pomts A, B, D, and so forth, on the airfoil surface may be 


First-family (Cj) 
Mach lines. 


Streamlines 


Mo 


FraunE 30.—Schematic diagram of shock-expansion method for cal- 


culating the flow field about an airfoil. 


calculated in the usual manner. If point B is chosen close 
to point A, the first-family (C,) Mach line connecting B to 
point C on the shock wave may be considered straight and 
inclined at an angle to the free stream equal to Вв--ӛв. 
Similarly, the segment of the shock wave AC may be con- 
sidered straight and inclined to the free-stream direction at 
an angle of c4. Thus, the physical coordinates of point C 
may be easily calculated. Since 5 is assumed constant along 
first-family Mach lines, бс is equal to ôs. All fluid properties 
at point C may be calculated from this known value of ёо 
with the oblique-shock-wave equations. In a similar 
manner, the segment DE of the next first-family Mach line 
is considered straight and inclined at Bp-+édp, and the 
streamline joining points C and E is considered straight and 
inclined at 8g. The physical coordinates of point E are 
therefore easily obtained. Since the flow along streamlines 
downstream of the shock wave is isentropic and since the 
pressure is also assumed constant along the first-family 
Mach lines (i. e., Ре--Рр), the fluid properties at point E 
are rendily obtained from the known properties at point C, 
using the isentropic flow relationships. The construction 
of the remainder of the flow field follows in a similar manner. 

As was discussed previously, the assumption that 6 is 
constant along first-family Mach lines is an additional con- 

15 It is clear, of course, that an average-value procedure could also be employed. In fact, 
for the shock-expansion method, an average-value procedure requires very little additional 
computation since the fluid properties for a system of points can be obtained independently 
of their physical coordinates. (The difference in the average- or initial-value procedyres 
appears only in computing the physical coordinates.) Thus, the slopes at both ends of a 


line segment are known before the line is added to the construction. No iteration, as re- 
quired іп characteristic solutions, 1s necessary in this case. 


dition which, in general, overdetermines the flow field. It is 
possible, therefore, to calculate two values of the shock-wave 
angle at each point on the shock, one assuming $ is constant 
and one assuming the pressure is constant. ‘These two values 
will differ slightly and it can be shown that they will bracket 
the correct value that would be given by the method of 
characteristics. It is also relatively easy to show that if 
the change in ô (ог p) along Cı given by the corresponding 
characteristic solution is small," the true value of the shock- 
wave angle lies just midway between the two values given by 
the different assumptions. It is apparent then that a 
closer approximation to the shock-wave shape can be easily 
obtained by simply averaging the values of shock-wave 
angle determined by assuming ё is constant and by assuming р 
is constant. The increase in accuracy is illustrated in figure 31 
for the biconvex airfoil at Mo= о. The averaging procedure 
gives a shock-wave shape that is closer to that given by the 
method of characteristics by 60 to 80 percent. The increase 
in computation time (at least for ideal-gas flows) is negligi- 
ble. In this regard, it is interesting to note that the shock- 
expansion solutions require less than 20 percent the computa- 
tion time of the characteristic solutions. 

The shock-expansion method is applicable to the determi- 
nation of the flow not only in the region adjacent to the 
airfoil (whether concave or convex), but also in the region 
downstream of the airfoil; hence, it may, for example, prove 
useful in downwash studies and the like. 
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FIGURE 31.—Shape of the shock wave for 10-percent-thick biconvex 
airfoil section at «—0? and Afj= о. 


п The necessary assumption here is that (Fe and (4), are nearly equal. This is 
. 3 


equivalent to the condition that pe із much less than one. (Bee eq. (40).) 
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APPENDIX С ` | 


_APPROXIMATE SOLUTION FOR PRANDTL-MEYER FLOW OF 
A CALORICALLY IMPERFECT GAS 


The following solution is obtained with an analysis similar 
to that used in Meyer’s original paper (ref. 23). А schematic 
didgram of the subject flow field is shown in figure 32. It is 
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Figures 32.—Schematic diagram of Prandtl-Meyer flow around a corner. 


evident that the change in flow-inclination angle for Prandtl- 
Meyer flow can be written as follows: 


óy — ó—(B8 n — B)-- (ew — o) 


Since the flow is isentropic, a given value of the local pressure 
will determine the Mach angle, 8. The problem, then, is to 
evaluate the angle, w. To this end, the velocity components 
tangential and normal to the first-family Mach lines may be 
expressed in the usual manner in terms of а potential o, thus, 


(C1) 


0 Р 5 
1 (C2) 
_ 1де 


It is clear, however, that these components are functions of 

w only; hence, it is convenient to define a new velocity poten- 

tial which is a function of w alone. Such a potential is 

(0)=5 (C4) 

The velocity components may then be written in terms of 
this new potential. 


u= (C5) 
a=, ` (C6) 

The resultant velocity is given by the expression 
10557 | (C7) 
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Equation (15), for conservation of energy, may be written 
in terms of the local temperature as follows: 


eT 


V2 (2%) RT+2RT (= i 5 (C8) 


The constant, A, is evaluated at the conditions existing up- 
stream of the expansion region; namely, 


24 O/T ny ا‎ 
Ata V+ ( 28) RT y+ 2RT» (аа) (09 


Equations (C7) &nd (C8) are then combined to yield 


^ p,— —2RT (253 فل‎ aq) +A (C10) 
Yi 1 eT —1 
It was shown previously, however, that 
v=yRT (C11) 


Equations (C6), (C10), and (C11) may therefore be combined 
to obtain the following relationship: 


2 0/T 
sso | 12-2 (2 tatu) (с> 
or 
үскі, 2 Гу Үс. ШЕЛ 9/T } = з 
с ен m rail [DC € ры ) атт = 
(C13) 
From the imperfect-gas relationship for we have 
ет 
те. 14-0 1) (T) тту (607—1)? (014) 


am 
+4 (т) e 


By substitution of this relation into equation (C13) there 18 


then obtained 
2 yi 1 2 8 Iz 2 
ал. LE) i P) 
where 
в/т 
т روبنز‎ RE | tot | 
"($- al : (016) 
T) gr—i LUI ال‎ 
yı (el —1y 
Now 
YT c 
п БЕСТІГІ (C17) 


For every value of Т/@ there is thus a particular valuo of 
q3/y,R6@. The function F(0/T) is therefore uniquely deter- 
mined for any value of а? since 7:26 is, of course, a constant. 
With this point in mind, let 


F(0/ T) — G(a? -y.R0) (C18) 
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a2/y,R8 


FIGURE 33.—Accuracy of approximation used in Appendix C to obtain 
solution for Prandtl-Meyer flow of a calorically imperfect gas. 


Figure 33 shows G(a?/y,R@) plotted as a function of а?/у,ЁӨ. 
This curve is approximated with the following simple 
relation: 

0.14 


2 = _а?_ 4 — 
Ala |y ,R0) — 0.38 29810-71 a] y, RO (C19) 
a? 
for 0.18 ES EL 
and 
G(a?|y,R0)—0 (C20) 


for 0 22<0. 18 


ШШ (C19) is also plotted in figure 33 to show the accu- 
racy of this approximation. Consider first the case when G 
is given by equation (C19) which is written in the form 


бары RO) ka سدم‎ = (C21) 
where, obviously, 
£—0.38/y,R0 
—0.71 (C22) 
n= — 0.14 (у.б) 


Equation (C21) is substituted into equation (C15) and, with 
equations (O18) and (C6), the following expression results: 


za Ф,44- CELERE) ترج‎ pa t A?=0 (C23) 


In n to ару this equation the lowing substitutions Р 


аге made: 
рі EE +25) 45 Hla (42 n) (029 
S£ ( At m 
sin? v= > ни (C25) 
8E وج‎ 
gin? كم‎ (C28) 
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and 


D^ cost +(r4) (cen 


سين 


Yi 
Equation (C23) then reduces to 


ب فى 3 


4 sin? + _ 4 sin’ y 


cos? TÔT 3ں‎ -- —— -p Ту 


4. сов у 
D 


Cost r(rw)*-+ =0 (C28) 


This equation is solved for 7,, thus, 


To =. ee (cos т —cos >)! 
cos rdr 
d= үе т 2 (cos т—сов ,) 2 


This expression is readily integrated to obtain the following 
equation relating w to the local velocity: 


(C29) 
or 


(C30) 


о-шн- Үр (2 [E(k,z) — E(k,z)] —[F(k,2)—F(k,2n)] } 


where 
Е eliptic integral of the second kind 
Е elliptic integral of the first kind 


k sin 5 (modulus) 


(C31) 


sm? 2) (amplitude) | 
The procedure for calculating corresponding values of the 

pressure, р, and the deflection angle, 5, із straightforward 
with the aid of the preceding equations and may be sum- 
marized as follows: 

1. Calculate A? (eq..(C9)) 

2. Calculate D? (eq. (C24)) : 

3. Calculate » (eq. (C25)) 

4. Assume a value of T, less than Ту 
5. Calculate р (eqs. (16) and (17)) 
6 
7 
8 


N 


sin~! 


. Calculate V? (eq. (С8)) 
. Calculate y (eq. (C14)) 
. Calculate а? (or Ф,2) (eq. (C11)) 
9. Calculate M and, in turn, 8 from V and a 

10. Calculate wu? (or &? (eq. (C7)) 

11. Calculate 7 (eq. (C26)) 

12. Calculate w (eq. (C31)) 

13. Caleulate 6 (eq. (C1)) 

This procedure is followed so long as the quantity а? [-y.R0 
is greater than 0.18. (This is equivalent to the temperature 
being greater than approximately 1,000? R.) For values 
of a?|y;R8 less than or equal to 0.18 (or temperatures less 
than about 1,000? R.), Gis set equal to zero (see eq. (C20)) 
In this case, equation (C15) reduces to the same form as for 
an ideal gas and, therefore, the well-known ideal-gas rela- 
tionships can be used. 


è 
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APPENDIX D 


EVALUATION OF (1/Ky)(dP/dW), K,/Kw, Y, AND к 


For an ideal gas, the Mach number, Mach angle, shock- 
wave angle, and pressure ratio can be calculated at the leading 
edge, using the standard Rankine-Hugoniot shock relations 
and utilizing the free-stream Mach number and leading- 
edge deflection angle. With these flow parameters known, 
the only terms to be determined іп equations (38), (39), 


(40), (43), and (44) are c4 and (% 4 2). These derivatives 


are easily obtained from the shock relations. (See ref. 9.) 
ару  (dP[de), 
(T). (айа). ir 
where 
_ 2y Му sin 2с Мо? sin 2с 
(4). ovi (94) 
апа 


4) ces ro 1) cos? (c —ô)+ sin 2(c —8) cot 2с (03) 


By again € the standard forms of the. Rankine- 
Hugoniot shock relations, it is possible to transform equations 
(38), (39), and (43) (given in the analysis) into 


E а 
mee сов: В +з) 


1 аР vtaent cos? t o? зір? о 
Ky dW sin? В | tan? 1 cos? By, l1 1 aa) 
tan? B 2 \cos? ¢ М? sin? о 
2y Mè sin? 3 | 
| +1 Van) 
for the surface pressure gradient, 
tan? с 
К, 4cost| tant 1 Cd ТҮН ЭБИН 1 as) 
tan? В ' 2 \cos?¢ Mo? sin? 
for the shock-wave curvature, and 
cos? B 1 
es 1+2 cos? ل‎ Weare) tan ¢ (D6) 
— | tan? P 1 (сов? В 1 


бал В 
ten! 72 cos E" M snis sin! ims.) 


for the surface-pressure-gradient ratio. ‘These equations are 
similar to those given by Schaefer in reference 8 and require 
Jess work to compute than equations (38), (39), and (43). 
For calorically-imperfect-gas flows, the standard shock 
relations are obviously not applicable. In this case, the flow 
parameters at the leading edge can be evaluated using the 
oblique-shock-wave equations previously presented. (See 
eqs. (14), (15), (18), (19), and (20).) Since the primary 
variable in these equations is temperature, the required 
362 


derivatives are most easily determined by employing the 


temperature as a parameter.’ ‘Thus, 
(%)- . (dP[d T). (D7) 
dé = (dsJdT), 
EM d (dé/dT) 
0 8 é 
2.) Guat (D8) 


Differentiation of equation (18) yields 
(2) Р p T(P— вала — М.Т, C pv oT) 
ат “ 


Т.+ТуР? 
(09) 
where from equation (15) 
dM MT, dy ‚М 
dT ub 274 2 t 
O/T. f Yt 
JT» It 
d (a i) + | OM) 
and from equation (14) 
ar ` 2 1) (y:—1) 575 вт, — 1 > 


The derivative (sr) can be evaluated from equation (20) 


dP 

yoM tan с t 

2 =sin’ | sirens (ар) ele 
e (P—1) 


зма pai) sec? с (Sr). | (D12) 


where from equation (19) 
(4) - 
— MAT, 
Yo 2 dM 1 
csc 2o (2 MT. 7) 06 jT dT TT 


(2) 1 


1 dy 
Ys dT 


tang ^5 ЖЕЗ) 


(D13) 
and from equation (8) 
(+4 2.) = Т. 2, 6 -»| (D14) 
The procedure for calculating = ФР, Ke, y, and xis straight- 


forward with the aid of the preceding equations and may be 
summarized as follows: 
1. For any Мо, To, бу, and X, choose a value of Т,/Ту. 


18 The forms of some of the derivatives presented havo been somewhat simplified from the 
forms originally presented in reference 17. 
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(Гір, 6 or the ideal-gas relations will provide an initial 
estimate.) 


. Calculate у, (eq. (14)) 

. Calculate М, (eq. (15)) 
. Calculate P, (eq. (18)) 
. Calculate рг/рг (eq. (8)) 
. Calculate с (eq. (19)) 

. Calculate 8, (eq. (20)) 


If this value of à, is not close enough to the desired value of 
dy, iterate, choosing a new Т,/ To. 


“IO Ov i C to 


8. Calculate 27, (sq. (D11) 

9. Calculate 2и (eq. (D10)) 
Calculate (27), (eq. (D9)) 
11. Caleulate (2% (eq. (014) 
12. Calculate (тт), (eq. (D13)) 
13. Calculate (Sr), (eq. (D12)) 
14. Calculate (е _(eq. (D7)) 


906/0 و‎ 
15. Calculate 55/00, (eq. (40)) 


1 dP 
16. Calculate кат (ед. (38)) 


17. Calculate y (eq. (43)) 
18. Caleulate (25) (eq. (D8) 


19. Calculate Же (eq. (39)) 
20. Calculate к (eq. (44)) | 
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TABLE I—FUNCTIONS FOR SLENDER-AIRFOIL METHOD 
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TABLE IL—SURFACE PRESSURE GRADIENT 
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TABLE HI.—SURFACE-PRESSURE-GRADIENT RATIO, y 
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TABLE IV.—SHOCK-WAVE CURVATURE, K,/K. 
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TABLE VI.—SURFACE PRESSURE GRADIENT, SURFACE- 
. PRESSURE-GRADIENT RATIO, SHOCK-WAVE CURVA- 
TURE, AND SHOCK-WAVE-CURVATURE RATIO FOR A 


CALORICALLY IMPERFECT, DIATOMIC GAS 
[To= 500° Rankine] 
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* Values of the temperature greater than 5,000? Б. downstream of the shock wave. 4 
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